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1.  EXECUTIVE  SUMMARY 


A  goal  of  the  Balanced  Technology  Initiative  (BTI)  on  Smart  Weapons  Operability 
Enhancement  (SWOE)  is  to  model  the  infrared  (IR)  radiant  field  from  complex  natural 
backgrounds.  The  degree  of  realism  and  model  detail  must  be  commensurate  with  the  spatial 
and  temporal  resolution  of  operational  and  developmental  sensors.  Thus,  the  background  model 
with  its  various  scene  elements  must  treat  the  basic  phenomena  affecting  the  radiation  signature 
to  the  required  resolution. 

The  overall  goal  of  this  Spectral  Sciences,  Inc.  (SSI)  effort  is  to  provide  a  preliminary 
study  that  quantifies  the  effect  of  cloud  shadows  on  surface  IR  radiance  and  thereby  helps 
determine  the  importance  of  including  shadowing  in  the  SWOE  model.  The  diurnal  variation  of 
surface  ground  temperature  has  been  modeled  for  an  "ideal"  dry  sandy  soil  with  12  cm 
vegetation  and  covered  by  various  stationary  cloud  fields.  Changing  cloud  shadows  result  from 
the  sun’s  movement  across  the  sky.  The  information  resulting  from  this  study  includes  (1) 
ground  temperature  response  to  shadowing  by  finite  clouds  and  (2)  diurnal  ground  temperature 
dependence  on  fractional  cloud  cover. 

The  general  approach  for  determining  ground  temperature  involves  three  steps: 

•  ATMOSPHERIC  RADIANCE:  MODTRANO >2)  with  its  multiple  scattering 

option  is  run  to  compute  the  diurnal  variation  of  downward  radiant  flux  and 
solar  illumination  for  cloud-free  and  semi-infinite  cloud  conditions; 

•  CLOUD  SHADOWING:  Cloud  fields  consisting  of  regularly  spaced  prolate 
ellipsoidal  clouds  are  constructed,  and  the  temporal  variation  of  solar 
illumination  is  determined;  and 

•  SURFACE  TEMPERATURE:  The  above  information  is  input  to  a 

time-dependent  energy  model  for  soil  terrains  with  a  simple  vegetation  layer 
to  determine  surface  temperature. 

Figure  I  illustrates  the  computed  variation  in  ground  temperature  for  40%  cover  by 
2.4-3.0  km  altitude  clouds.  The  model  is  run  for  at  least  one  full  day  to  reduce  the  bias 
created  by  initialization  of  the  soil  temperature  profile.  As  a  general  trend,  the  ground  warms 
after  sunrise  (*4:30  a.m.  on  June  25  at  Hanscom  AFB)  and  cools  during  late  afternoon  and 
into  the  night.  The  interesting  features  in  this  curve  are  the  sharp  rises  and  drops  which 
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TIME  (HRS  RFTER  MIDNIGHT) 

Figure  I.  Diurnal  Variation  of  Ground  Temperature  for  70%  Cover  by  Finite  Clouds 
Between  2.4  and  3.0  km  Altitude. 

correspond  to  breaks  in  the  cloud  cover.  The  discontinuities  in  solar  irradiance  produce 
changes  of  up  to  2K  in  one  minute,  and,  within  ten  minutes,  the  ground  temperature  changes 
by  more  than  4K. 

The  model  also  predicts  the  different  ground  temperature  histories  that  arise  from  varying 
cloud  cover  while  keeping  all  other  model  inputs  constant.  Curves  of  ground  temperature  for 
10,  40,  70,  and  100%  cover  by  2.4-3.0  km  clouds  were  computed  relative  to  the  cloud-free 
ground  temperatures;  the  results  for  a  24  hour  period  are  shown  in  Figure  2.  For  dry,  sandy 
soil,  full  cloud  cover  reduces  ground  temperature  by  more  than  10°C  at  midday  and  increases 
the  night  temperature  in  excess  of  1°C.  For  fractional  cloud  covers,  the  early  morning  ground 
temperatures  vary  smoothly  between  the  clear  sky  and  complete  cloud  cover  limits.  This  is  not 

the  case  during  the  day.  The  ground  temperature  responds  quickly  to  the  abrupt  changes  in 

solar  illumination.  Figure  2  also  illustrates  that  the  limiting,  stable  ground  temperatures  are 
approached  quite  slowly.  Consider,  for  example,  the  40%  cloud  cover  curve  between  36.3  and 

39.3  hours  after  midnight.  During  this  entire  period,  the  ground  is  solar  illuminated  and  the 
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Figure  2.  Ground  Temperature  Differences  Between  Cloud  and  No  Cloud  Scenarios. 

The  Cloud  Field  is  Constructed  from  Regularly  Spaced  Prolate  Ellipsoidal 
Clouds  Between  2.4  and  3.0  km  Altitude. 

skyshine  is  essentially  the  same  as  the  cloud-free  case.  Even  so,  over  two  hours  pass  before 
the  ground  temperature  is  stabilized.  Thus,  the  model  indicates  that  although  the  surface  reacts 
quickly  to  changes  in  solar  illumination,  the  sub-surface  layers  react  quite  slowly. 

The  major  results  of  this  work  are  a  quantitative  verification  that  cloud  shadowing 
significantly  affects  the  radiant  field  from  complex  natural  backgrounds  and  that  predictive 
modeling  can  describe  the  temporal  and  spatial  responses  of  the  soil  surface  temperature.  Thus, 
this  study  demonstrates  the  necessity  of  including  cloud  shadowing  in  the  SWOE  scene  model. 
This  effort  has  also  produced  a  series  of  modules  which  can  serve  as  a  basis  for  development 
of  a  comprehensive  cloud  shadowing  model.  Model  upgrades  and  validation  will  be  needed  to 
insure  that  SWOE  performance  requirements  are  met. 
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2.  CLOUD  SHADOWING  MODEL 


In  order  to  quantify  the  importance  of  cloud  shadowing  for  the  SWOE  natural 
background  model,  SSI  has  developed  a  preliminary  cloud  shadowing  model.  The  model  consists 
of  three  modules: 

•  DIURNAL  -  Calculates  the  diurnal  variation  of  radiation  flux  incident  on  the 
ground, 

•  CLOUDS  -  Constructs  cloud  fields  and  determines  direct  solar  transmittance, 
and 

•  TERRAIN  -  Solves  a  time-dependent  energy  budget  equation  for  surface 
temperature. 

The  following  subsections  discuss  the  details  of  these  models,  their  implementation,  and  results. 

Surface  temperature  predictions  due  to  cloud  shadows  are  naturally  very  sensitive  to  input 
physical  conditions.  Table  I  is  a  complete  list  of  model  parameters.  The  model  includes 
diurnal  variations  in  solar  irradiance  and  scattering,  and  calculations  were  performed  for  differing 
cloud  conditions.  Since  this  preliminary  study  is  to  demonstrate  quantitatively  the  effect  of 
changing  cloud  cover  on  the  radiation  signature,  representative  soil  and  vegetation  parameters 
were  chosen  to  correspond  to  a  scenario  in  which  the  temperature  response  is  expected  to  be 
accentuated.  The  vegetation  parameters  correspond  to  summer  barley,  and  surface  conditions 
reflect  a  fairly  dry  and  sandy  soil. 

2.1  Diurnal 

DIURNAL  is  a  driver  routine  which  repeatedly  calls  M0DTRAN(L2)  to  determine  the 
diurnal  variation  in  solar  irradiance  and  downward  flux  of  skyshine  radiation  on  the  ground. 
Inputs  include  a  reference  MODTRAN  input  stream  along  with  time  step  information. 
Beginning  at  sunrise,  DIURNAL  executes  MODTRAN  at  each  time  step  with  solar  direr  ion 
and  flux  determined  by  time  of  day,  day  of  year,  longitude  and  latitude. 

The  only  difference  between  the  MODTRAN  input  streams,  other  than  time  of  day,  is  an 
adjustment  made  to  the  reference  temperature  and  water  vapor  profiles  in  the  boundary  layer 
(<0.L5  km)  to  insure  consistency  with  surface  values.  Following  Barnett,^)  diurnal  variations  in 
surface  (1.5  m  above  ground)  temperature,  T,  and  relative  humidity,  R,  are  parameterized  in  the 
cloud  shadowing  model,  with  extreme  values  occurring  at  sunrise,  ts,  and  2:00  p.m.: 


TABLE  I.  CLOUD  SHADOWING  MODEL  INPUT  PARAMETERS 


1  EXTRA-TERRESTRIAL  SOLAR  IRRADIANCE 

|  pjm;  Calculated  by  MODTRAN  for  day  176  (June  25) 

|  Direction  Calculated  every  20  minutes  on  day  176 

j  at  42°28’N,  71°1TW  (MODTRAN) 


ATMOSPHERE 

T,  P  &  Molecular  Profiles  US  STANDARD 

Aerosol  Profile  23  km  Visibility  Rural  Extinction 

Model  (MODTRAN) 


CLOUD  FIELD 

Distribution 

j  Absorption 
j  Scattering 


Regularly  Spaced  Prolate  Ellipsoidal 
Clouds  at  0.7  and  2.7  km  Altitude 
MODTRAN  STRATUS  and  ALTOSTRATUS  Models  j 
MODTRAN  STRATUS  and  ALTOSTRATUS  Models  j 


BOUNDARY  LAYER  (1.5  m 

Wind  Speed 
Air  Temperature 
Relative  Humidity 


above  ground  averages) 

2.0  m/s 

Varies  Smoothly  between  278  and  300  K 
Varies  Smoothly  between  40  and  45% 


j  VEGETATION 

!  Fractional  Cover 
Long  Wave  Emissivity 
Short  Wave  Albedo 
Height 

SOIL 

Initial  Temperatures 

Conductivity 

Diffusivity 

Saturation 

Heat  Transfer  Coefficient 
Long  Wave  Emissivity 
Short  Wave  Albedo 


0.3 
0.95 
0.20 
0.12  m 


Evenly  Spaced  between  Surface  Air  Temperature 

and  290K  at  1  m  Depth 

1.1  W  or1  K'1 

0.004  cm^/s 

25% 

0.012 

0.95 

035 


T  = 

To 

2  s  t  s  ts 

(1) 

T  = 

T0  +  AT  sin(^ 

ts  <  t  <  14 

(2) 

T  = 

(T0  +  aT)^26'1)/12  Tq^'14)/12 

14  s  t  s  26 

(3) 

R  = 

R0 

2  s  t  s  ts 

(4) 
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R 


ts  <  t  s  28-ts 


(5) 


=  R0  -  AR  sin(y 


-ts 

14-t 


R  =  R0 


28-ts  <  t  <  26 


(6) 


where  t  is  the  local  time  in  hours  after  midnight  (2  to  26),  T0  and  R0  are  the  temperature 
and  relative  humidity,  respectively,  at  sunrise,  AT  and  AR  <.re  the  diurnal  range  of  temperature 
and  relative  humidity,  respectively,  and  sunrise  is  assumed  to  occur  a  minimum  of  two  hours 
after  midnight. 

The  five  separate  MODTRAN  runs  which  are  performed  at  each  time  step  are: 

•  Clear  Sky  Direct  Solar  Irradiance, 

•  Clear  Sky  Short-Wave  (0.3-3.3  pm)  Skyshine, 

>  Clear  Sky  Long-Wave  (3.3-200.0  pm)  Skyshine, 

•  Clouded  Sky  Short-Wave  (0.3-3.3  pm)  Skyshine,  and 

•  Clouded  Sky  Long-Wave  (3.3-200.0  pm)  Skyshine. 


MODTRAN  has  been  modified  to  output  skyshine  (downward  radiant  flux)  which  is  calculated 
as  part  of  the  multiple  scattering  routine.  Separate  calculations  are  performed  in  the  long-  and 
short-wave  to  accommodate  different  surface  albedos.  The  DIURNAL  assumption  that  the 
surface  albedo  is  independent  of  both  angle  and  wavelength  within  each  spectral  region  is 
consistent  with  the  present  version  of  the  energy  budget  model.  Cloud  calculations  for 
determining  skyshine  assume  100%  cloud  cover.  Since  skyshine  is  a  fairly  diffuse  source,  the 
cloud  shadowing  module  interpolates  the  downward  radiant  flux  between  the  clear  sky  and  100% 
cloud  cover  limits  for  input  to  the  energy  budget. 

The  DIURNAL  calculations  were  performed  for  an  observer  located  at  Hanscom  AFB 
(42°28’N,  71°17W),  Table  I.  Radiation  flux  was  determined  each  twenty  minutes  for  a  24  hour 
period  starting  from  sunrise  (=*4:30  a.m.)  on  day  176  (June  25).  The  results  were  interpolated 
temporally  for  input  into  the  heat  flow  equation  (Subsection  2.3).  The  1976  US  Standard 
atmosphere  and  the  MODTRAN  23  km  visibility  rural  extinction  aerosol  model  defined 
atmospheric  conditions.  Predictions  for  two  cloud  types,  low  altitude  stratus  (0.33-1.00  km)  and 
mid  altitude  altostratus  (2.4-3.0  km),(4)  were  calculated. 
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Surface  irradiances  from  the  DIURNAL  calculations  are  illustrated  in  Figure  3.  A 
number  of  observations  are  noteworthy: 

#  •  At  the  maximum  in  the  cloud-free  sky  solar  irradiance  curve,  corresponding  to 

a  19°  zenith  angle,  almost  60%  of  the  incident  flux  is  transmitted  directly. 

With  complete  cloud  cover,  essentially  no  unscattered  solar  radiation  reaches 

.  the  ground. 

•  At  night,  the  downward  radiation  flux  is  enhanced  by  both  the  emission  of 

thermal  energy  by  clouds  and  the  clouds’  reflection  of  earth’s  thermal 

emission.  For  warm,  low  altitude  clouds,  this  effect  is  greatest. 

•  Solar  scatter  contributes  to  skyshine  during  the  day.  Thin  clouds  tend  to 

increase  the  amount  of  solar  radiation  scattered  downward,  but  as  the  cloud 

optical  depth  becomes  large,  most  solar  radiation  is  either  reflected  upward  or 

absorbed.  Since  MODTRAN  extinction  coefficients  for  stratus  clouds  are 
approximately  60%  of  those  for  altostratus,  the  stratus  clouds  lead  to  a 
greater  contribution  from  solar  scatter. (5) 

•  During  early  morning  and  late  afternoon,  cloud  depth  along  the  solar  path 

becomes  quite  large,  and  the  solar  scatter  contribution  to  skyshine  drops.  In 
fact,  the  total  skyshine  for  altostratus  clouds  is  below  the  clear  sky  value 

between  5:50  and  9:30  a.m.,  and  again  between  2:30  and  6:10  p.m. 


TIME  (HRS  RFTER  MIDNIGHT) 

Figure  3.  Diurnal  Variation  of  Surface  Irradiance  Components  Calculated  by  DIURNAL. 

Skyshine  Includes  Scattered  Solar.  The  Dashed  Curves  Include  Complete 
Cover  by  Stratus  (0.33-1.00  km  altitude)  and  Altostratus  (2.4-3.0  km  Altitude) 
Clouds.  The  Predictions  are  for  an  Observer  at  Hanscom  AFB  on  Day  176. 
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•  At  midday,  the  total  skyshine  contribution  from  the  clouds  at  0.7  km  altitude 
actually  exceeds  the  clear  sky  direct  solar  component.  Of  course,  this  does 

not  mean  the  clouded  sky  is  warmer  than  the  cloud-free  sky,  the  total 

downward  flux  in  the  cloud-free  sky  scenario  is  the  sum  of  the  solar 

irradiance  (776  W/m^)  and  skyshine  (54 8  W/m^)  contributions,  which  are 

considerably  more  than  the  stratus  cloud  skyshine  (796  W/m^). 

2.2  Clouds 

The  routine  CLOUDS  determines  the  diurnal  variation  in  transmittance  of  direct  solar 
irradiance  through  a  stationary  cloud  field.  For  this  study,  cloud  fields  were  defined  as  a 
lattice  of  prolate  ellipsoidal  clouds,  each  centered  at  a  fixed  altitude  above  the  curved  earth  and 

oriented  so  that  the  local  vertical  is  parallel  to  the  major  semi-axis.  Figure  4  is  a  polar  plot 

illustration  of  the  cloud  field  as  seen  by  a  ground  observer.  In  this  cloud  lattice,  the  observer 
always  has  a  cloud  directly  overhead.  The  figure  was  generated  for  40%  cover  by  clouds  of 
width,  base  and  top  equal  to  1.0,  2.4,  and  3.0  km,  respectively.  (Throughout  this  report,  cloud 

cover  is  defined  as  the  percentage  of  the  earth’s  surface  for  which  clouds  are  directly  overhead; 

this  is  different  from  the  percent  of  blue  sky  blocked  by  clouds  from  a  ground  observer’s 
perspective.)  Figure  4  is  equivalent  to  an  all-sky  image;  the  radial  distance  from  the  center  is 
proportional  to  zenith  angle  with  the  edge  corresponding  to  the  horizon,  and  the  azimuthal 
angle  is  defined  relative  to  the  compass  directions.  The  periodic  cloud  structure  was  selected  to 
lessen  the  severity  of  the  assumption  that  skyshine  can  be  interpolated  between  the  cloud-free 
and  100%  cloud  cover  limits. 

Additional  simplifying  assumptions  were  instituted  in  this  preliminary  cloud  shadowing 
model.  The  sun  is  treated  as  a  point  source,  and  the  clouds  are  homogeneous.  Initially, 
optical  properties  for  the  finite  cloud  were  defined  with  the  MODTRAN  model,  and 

transmitt  an  ces  were  calculated  with  the  SSI  cloud  model.  The  CLOUDS  module  calculates  the 

diurnal  variation  of  geometric  cloud  depth  along  the  solar  path,  as  shown  in  Figures  5  and  6. 

Since  cloud  extinction  coefficients  are  so  large  in  the  short  wave  (>  88  km'1  for  altostratus  and 

>  53  km'1  for  stratus(4)),  the  transmittance  rapidly  falls  to  zero  whenever  clouds  intercept  the 
solar  path.  Therefore,  the  spectral  integral  over  the  product  of  solar  irradiance  and  cloud 
transmittance  was  eliminated,  and  the  solar  irradiance  set  to  zero  whenever  clouds  intersect  the 
solar  path  to  the  ground. 
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NORTH 


SOUTH 

Figure  4.  Schematic  Illustrating  the  Stationary  Cloud  Field.  The  Polar  Plot  Solar 

Trajectory  is  for  a  Ground  Observer  Located  42°28’N,  71°17’W  on  Day  176. 
Distance  from  Center  is  Proportional  to  Zenith  Angle  with  the  Edge 
Representing  the  Horizon.  Cloud  Cover  is  40%  and  the  Clouds  are  Prolate 
Ellipsoids  with  Base  =  2.4  km,  Top  =  3.0  km,  and  Width  =  1.0  km.  Horn- 
Glasses  are  Drawn  Each  Hour  Past  Sunrise  with  Sunrise  Beginning  in  the 

East  at  4:30  am.  Local  Time. 

A  lattice  structure  for  the  cloud  fields  was  chosen  in  order  to  be  more  consistent  with 
the  semi-infinite  clouds  used  to  calculate  the  skyshine  irradiance.  The  MODTRAN  altostratus 
cloud  model  has  a  base  of  2.4  km  and  top  of  3.0  km;  the  values  for  the  low  altitude  stratus 

are  0.33  km  and  1.00  km,  respectively.  All  clouds  were  given  a  width,  W,  of  one  kilometer. 

The  distance  between  nearest  neighbor  clouds,  D,  was  calculated  from  the  equation 

D  =  -y  iTn/F  ,  (7) 
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Figure  5.  Diurnal  Variation  of  Geometric  Cloud  Depth  Along  the  Solar  Path  for  Three 
Fractional  Covers  of  Clouds  at  2.4  km.  The  Cloud  Scene  Consists  of  a 

Square  Lattice  of  Prolate  Ellipsoidal  Clouds  Each  with  Base,  Top,  and  Width 
Equal  to  2.4,  3.0,  and  1.0  km.  Respectively.  The  Ground  Observer  is  located 
42°28’N,  71°1TW,  and  Time  is  Measured  from  the  Beginning  of  Day  175 
(4:00  a.m.  to  8:00  p.m.  on  Day  176  at  Hanscom  APB). 
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Figure  6.  Diurnal  Variation  of  Geometric  Cloud  Depth  Along  the  Solar  Path  for  Three 
Fractional  Covers  of  Clouds  at  033  km  and  for  the  Same  Conditions  as 
Figure  5. 
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where  F  is  the  fraction  of  the  earth  covered  by  the  cloud  field.  Note  that  F  is  smaller  than 
the  fraction  of  sky  containing  clouds  from  the  perspective  of  a  ground  observer.  Calculations 
were  performed  for  F  values  of  10,  40  and  70%. 

Cloud  shadowing  was  modeled  for  day  176,  June  25.  The  sun’s  trajectory  as  seen  from 
the  ground  at  Hanscom  AFB  is  shown  in  Figure  4.  Note  that  the  F=0.4  curve  in  Figure  5  is 
indeed  zero  whenever  the  solar  path  passes  between  clouds  in  Figure  4.  The  cloud  directly 
overhead  always  occults  the  sun  around  noon  for  the  low  altitude  cloud  lattice,  so  the  three 
factional  cloud  covers  in  Figure  6  are  the  same  at  36  hours. 

2.3  Terrain 


Calculations  of  time-dependent  soil/vegetation  temperature  were  performed  by  the 
TERRAIN  module.  This  module  is  based  on  Deardorffs  terrain  model(^)  with  the  single 
additional  limitation  that  the  soil  moisture  content  is  constant.  TERRAIN  uses  a  finite 
differences  approach  to  solve  the  one-dimensional  heat  flow  equation  below  the  surface: 


3T  a^T 

at  =  K  dz2 

where 

K  =  soil  thermal  diffusivity, 

z  =  depth  below  the  surface, 

t  =  time,  and 

T  =  temperature  at  depth  z  and  time  t. 


(8) 


The  boundary  conditions  include  (1)  an  initial  (early  morning)  temperature  profile,  (2)  a 
constant  lower  boundary  temperature  (one  meter  below  the  soil  surface)  and  (3)  a  temporal 
surface  temperature  which  is  determined  along  with  the  vegetation  temperature  from  a  pair  of 
energy  budget  equations: 


Ss  + 

Rs  +  Hs  +  Es  +  Gs  =  0 

(9) 

Sy  + 

Ry  +  Hy  +  Ey  =  0 

(10) 
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where 


S  =  Short  wave  (solar)  illumination, 

R  =  Net  Long  wave  (thermal)  radiation, 

H  =  Sensible  heat  exchange, 

E  =  Evapotranspiration  heat  exchange,  and 

G  =  Conductive  transfer  to  the  soil. 

Subscripts  s  and  v  refer  to  soil  and  vegetation,  respectively. 

Beginning  with  an  initial  temperature  profile  for  the  501  homogeneous  soil  layer 
boundaries,  TERRAIN  increments  in  time,  solves  the  heat  flow  equation  to  determine  the 

temperature  profile  below  the  surface  (every  3  seconds),  and  then  uses  the  energy  budget 
equations  to  determine  the  soil  and  vegetation  temperature  (every  15  seconds).  The  model  is 
run  for  2  days  with  results  for  the  second  day  reported;  the  first  day  allows  the  soil 

temperature  profile  to  converge  to  its  periodic  limit.  Inputs  to  the  model  include  the  time 

dependent  skyshine  and  solar  irradiance  from  the  DIURNAL  and  CLOUDS  modules  along  with 
the  boundary  layer,  vegetation  and  soil  parameters  listed  in  Table  I. 

To  maximize  the  temperature  response,  a  fairly  dry,  sandy  soil  was  selected.  Low 
moisture  content  and  air  pockets  lower  the  overall  conductivity,  and  the  low  moisture  content 
also  decreases  evapotranspiration.  The  net  effect  is  that  absorbed  radiation  has  no  quick  outlet 
and  therefore  produces  rapid  temperature  changes.  The  soil  parameters  in  Table  I  are  based 
on  averages  listed  in  the  first  two  chapters  of  Boundary  Laver  Climates/6)  Soil  moisture, 

although  low  (25%  of  saturation),  is  in  the  acceptable  range  for  plant  growth.(6)  The  bare  soil 
heat  transfer  coefficient,  c,  is  calculated  from 

[ln(z/Zo)]2  (U) 

where  k  is  the  von  Karmon  constant  (0.4),  z  is  the  altitude  of  the  surface  wind  speed 
measurement  (1-5  m)  and  Zq  is  the  roughness  length  (0.0016  m  used  here;  compare  Table  2.2 
of  Reference  6). 

$  DeardorfPs  model(^)  allows  for  introduction  of  a  simple  vegetation  layer.  Vegetation 

height  and  fractional  cover  have  been  assigned  values  appropriate  for  summer  barley.^’2)  The 
vegetation  fractional  cover  is  defined  as  the  area  average  shielding  factor  associated  with  the 
degree  to  which  foliage  prevents  short-wave  radiation  from  reaching  the  ground.  Estimates  for 
short  wave  albedo  and  long  wave  emissivity  are  based  on  Table  1.1  of  Reference  6. 
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3.  CALCULATIONS 


p 


Calculations  were  performed  for  both  altostratus  (2.4-3.0  km  altitude)  and  stratus  clouds 

(0.33-1.00  km  altitude)  with  fractional  cloud  covers  of  0,  10,  40,  70  and  100%.  Results  are 

discussed  the  following  subsections. 

3.1  Surface  Energy  Budget 

Solutions  to  the  energy  budget  equations  show  the  response  of  vegetation  and  soil  to 

changes  in  incident  radiation.  In  this  subsection,  four  scenarios  are  considered:  foliage  and  soil 
energy  budgets  with  no  clouds  and  with  40%  cover  by  the  stationary  altostratus  cloud  field. 

The  diurnal  variation  of  the  foliage  energy  budget  is  driven  by  the  incident  radiation, 
which  consists  of  solar  flux,  skyshine  and  soil  thermal  emission.  Since  the  coverage  for  summer 
barley  is  only  30%, (5.7)  most  of  the  incident  flux  passes  through  the  vegetation  without  being 
absorbed.  For  clear  sky  conditions.  Figure  7  reveals  that  the  nighttime  incident  radiation  is 

more  than  balanced  by  the  vegetation’s  thermal  emission.  During  the  day,  both  solar 
illumination  and  convection  from  the  warm  air  add  energy  to  the  barley.  This  energy  raises 
the  foliage  temperature  and  causes  foliage  moisture  to  evaporate. 

Cloud  shadowing  produces  sharp  changes  in  the  solar  irradiance.  Figure  8  illustrates  the 
energy  budget  for  40%  cover  by  altostratus  clouds.  Each  valley  in  the  absorption  curve 
corresponds  to  a  cloud  shadow.  To  compensate  for  the  decrease  in  incident  radiation,  foliage 
thermal  emission  and  water  evaporation  decrease  and  sensible  heat  transfer,  i.e.,  convection  from 
the  surrounding  air  mass,  increases.  The  differences  between  nighttime  absorption  in  Figures  7 
and  8  is  less  than  the  differences  in  Figure  3  because  cloud  cover  is  only  40%  and  the  ground 
ts  warmer  under  cloud  cover  (i.e.  greater  soil  thermal  emission). 

Figure  9  illustrates  the  soil  energy  budget  for  this  cloud-free  scenario.  An  additional  flux 
term,  energy  conducted  through  the  soil,  is  present.  At  night,  the  top  soil  is  cooler  than  the 
lower  layers,  and  upwards  conduction  plus  the  downward  radiant  flux  keeps  the  surface  warm. 
During  the  day,  solar  irradiance  heats  the  surface,  causing  the  temperature  to  rise  and  reversing 
the  direction  of  the  soil  heat  flow.  Evapotranspiration  and  convection  are  of  minor  importance 
for  the  dry,  sandy  soil  under  cloud-free  conditions. 

Changes  in  the  soil  energy  balance  resulting  from  cloud  shadowing  are  shown  in 
Figure  10.  Consider  the  model  predictions  between  35  and  37  hours.  Within  seconds  of  solar 
occultation,  the  surface  temperature  begins  to  drop,  causing  a  sharp  soil  temperature  gradient 
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Figure  7.  Diurnal  Variation  of  the  Cloud-Free  Energy  Budget  for  Barley.  Absorption 
Includes  Solar  Irradiance,  Skyshine,  and  Soil  Thermal  Emission. 
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Figure  8.  Diurnal  Variation  of  the  Barley  Energy  Budget  Under  40%  Cover  by 
2.4-3.0  km  Altostratus  Clouds.  Absorption  Includes  Solar  Irradiance,  Skyshine, 
and  Soil  Thermal  Emission. 
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Figure  9.  Diurnal  Variation  of  the  Cloud-Free  Energy  Budget  for  Soil  Containing  Barley. 

Absorption  Includes  Solar  Irradiance,  Skyshine,  and  Vegetation  Thermal 
Emission. 
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Figure  10.  Diurnal  Variation  of  the  Energy  Budget  for  Soil  Containing  Barley  and 
Under  40%  Cover  by  2.4-3.0  km  Altostratus  Clouds. 
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and  reversing  the  direction  of  heat  conduction.  As  shadowing  persists,  the  temperature 
continues  to  drop,  but  the  magnitude  of  the  temperature  gradient  lessens  as  the  inner  soil 
layers  cool.  Energy  balance  is  maintained  by  reducing  thermal  emission,  evapotranspiration  and 
sensible  heat  transfer.  When  the  sun  reappears,  essentially  the  reverse  process  occurs;  a  sharp 
temperature  gradient  is  followed  by  a  reduction  in  its  magnitude  and  by  initiation  of  response 
from  the  other  energy  components. 


3.2  Diurnal  Surface  Temperatures 


The  long-wave  upward  radiant  flux,  I*,  from  soil  with  a  simple  vegetation  layer  has  the 

form 

I*  =  (l-ov)  [es  o  T*  +  (l-«s)  IA]  +  sv  [ev  o  T*  +  (l-ev)  IA]  (12) 

where 

av  =  Fractional  Cover  of  Vegetation, 
es  =  Soil  Surface  Emissivity, 

«v  =  Vegetation  Emissivity, 

Ts  =  Soil  Surface  Temperature, 

Tv  =  Vegetation  Temperature, 

a  =  Stefan- Boltzmann  Constant,  and 
IA  =  Long-Wave  Downward  Radiant  Flux. 


The  TERRAIN  module  calculates  soil  surface  and  vegetation  temperatures  so  that  scene 
radiances  can  be  determined. 


3.2.1  Cloud-Free  Conditions 


The  TERRAIN  model  predictions  for  the  diurnal  variation  of  vegetation,  soil  surface  and 
ground  temperatures  under  cloud-free  conditions  are  shown  in  Figure  11.  The  ground 
temperature  has  been  defined  as  a  weighted  average  of  the  soil  surface  and  vegetation 
temperature 


-  16  - 


(l-°v)  Ts  +  Oy  Tv 


(13) 


Also  shown  in  Figure  11  is  a  curve  of  the  surface  air  temperature,  which  is  input  to  the  model 
and  held  constant  for  these  scenarios. 

For  the  cloud-free  scenario,  the  barley  temperature  is  cooler  than  the  surface  air  at  all 
times  with  the  maximum  difference  approaching  3K  at  noon.  The  barley  remains  relatively  cool 
because  absorption  of  incident  daytime  solar  radiation  is  balanced  by  vaporization  of  the 
vegetation’s  moisture,  Figure  7.  The  soil  temperature  at  sunrise  is  about  IK  above  air 
temperature.  As  the  sun  begins  to  rise,  the  soil  heats  up  more  slowly  than  the  air.  Much  of 
the  initial  solar  irradiance  acts  to  warm  the  inner  soil  layers.  As  solar  irradiance  continues  to 
increase,  the  soil  surface  temperature  rises  more  sharply,  peaking  between  noon,  the  maximum 
in  solar  radiation  and  2  p.m.,  the  maximum  of  the  air  temperature.  Cooling  is  most  rapid 
between  4  and  7  p.m.,  almost  4K/hour  and  then  slows  as  the  conduction  of  heat  from  the 
inner  soil  layers  becomes  influential. 

3.2.2  Full  Cloud  Cover 

The  major  effects  of  complete  cover  by  2.4-3.0  km  altostratus  clouds  are  to  lower  the 

incident  radiation  during  the  day  by  occulting  the  solar  irradiance  and  to  increase  the  nighttime 

skyshine  by  forming  a  layer  of  insulation.  The  barley  responds  with  slightly  lower  temperatures 
during  the  day  and  slightly  warmer  temperatures  at  night,  Figure  12.  The  influence  on  the 

vegetation  temperature  is  not  large  because  the  solar  radiation  primarily  functions  to  lessen 
moisture  content  and  the  additional  nighttime  thermal  contributions  are  balanced  by  a  reduction 
in  convective  heating.  The  ..oil's  response  is  much  more  drastic.  The  maximum  soil  surface 
temperature  with  complete  cover  by  altostratus  is  15K  less  than  the  cloud-free  value.  From  2 
p.m.  to  2  a.m.,  the  soil  temperature  drops  at  a  relatively  slow,  almost  constant  rate  of 

LK/hour.  At  sunrise,  the  soil  temperature  with  full  cloud  cover  is  more  than  IK  warmer  than 
the  cloud-free  soil  temperature  -  a  reversal  of  over  16K  between  the  cloud-ffee  and  full  cloud 
cover  scenarios. 

3.2.3  Partly  Cloudy  Conditions 

The  response  of  the  soil  surface  and  foliage  temperatures  under  partly  cloudy  conditions 
with  2.4  -  3.0  km  altostratus  is  illustrated  in  Figures  13-15.  In  these  curves,  the  general  trend 
is  the  same.  As  soon  as  the  sun  is  occulted,  temperatures  drop  drastically,  and  when  the  sun 
reappears,  temperatures  abruptly  rise.  For  vegetation,  the  model  predicts  that  the  full  change 
occurs  instantaneously,  which  illustrates  a  limitation  of  the  model  in  its  present  form.  Since 
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Figure  11.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  Cloud-Free  Conditions. 
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Figure  12.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  100%  Cover  by  Altostratus  Clouds,  2.4-3.0  km. 
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Figure  13.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  10%  Cover  by  Altostratus  Clouds,  2.4-3.0  km. 
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Figure  14.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  40%  Cover  by  Altostratus  Clouds. 
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Figure  15.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  70%  Cover  by  Altostratus  Clouds. 

foliage  is  treated  as  a  single  layer,  there  is  no  mechanism  for  allowing  the  heat  to  properly 
dissipate  downward  through  the  vegetation.  For  soil  surface  temperatures,  there  is  also  a  quick 
response,  but  the  temperature  continues  to  change  for  over  an  hour.  Figure  16  shows  an 
enlargement  of  the  40%  cloud  cover  ground  temperature  curve  between  35.7  and  36.7  hours. 

The  temperature  changes  by  more  than  2K  in  the  first  minute,  and  by  a  full  5K  in  30  minutes. 

In  Figure  17,  ground  temperatures  for  10,  40,  70  and  100%  cover  by  altostratus  clouds 
(2. 4-3.0  km)  are  plotted  relative  to  the  cloud-free  temperatures.  At  night,  cloud  cover  provides 
an  insulation  for  the  atmosphere,  and  by  sunrise  ground  temperature  varies  smoothly  between 
the  clear  sky  and  semi-infinite  cloud  limits.  As  seen  in  Figure  3,  the  daytime  skyshine 
irradiance  from  cloud-free  and  100%  cloud  cover  by  2.7  km  altostratus  are  almost  the  same. 

Since  the  skyshine  under  these  conditions  is  relatively  independent  of  cloud  cover,  daytime 
temperature  variations  are  driven  by  the  direct  solar  irradiance,  and  the  ground  temperatures 
oscillate  between  the  cloud  cover  limits.  Figure  17  illustrates  the  considerable  lag  time  before  ■* 

the  soil  reaches  asymptotic  temperatures.  After  hours  of  solar  illumination,  previously  shaded 
ground  still  has  not  heated  to  the  clear  sky  limit.  Similarly,  regions  shaded  for  over  an  hour 
have  not  cooled  to  the  semi-infinite  cloud  limit. 
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Figure  16.  Enlargement  of  the  Ground  Temperature  Curve  from  Figure  14  During  the 
Passing  of  a  Finite  Cloud. 
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Figure  17.  Differences  in  Ground  Temperatures  Between  Altostratus  (2.4-3 .0  km)  Clouds 
of  Various  Fractional  Cover  and  Cloud-Free  Conditions. 
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3.2.4  Low-Altitude  Cloudy  Conditions 


An  additional  set  of  calculations  was  performed  for  low-altitude  stratus  (0.33-1.00  km) 
clouds.  Curves  of  the  diurnal  variation  in  soil  surface  and  summer  barley  temperature  are 
plotted  in  Figures  18-21.  Two  major  differences  exist  between  the  altostratus  and  stratus  clouds: 

•  For  a  given  fractional  cloud  cover,  the  lower  altitude  cloud  field  blocks  solar 
illumination  more  hours  of  the  day  (compare  Figures  5  and  6);  and 

•  Skyshine  from  warm  stratus  clouds  is  more  intense  than  that  from  altostratus 
(Figure  3). 

The  first  point  results  from  the  assumed  cloud  lattice  structure.  The  net  effect  of  these 
differences  is  that  the  ground  under  stratus  cloud  cover  is  warmer  than  under  altostratus  except 
when  solar  irTadiance  has  been  blocked  for  an  extended  period  of  time  by  the  low  clouds. 
Sudden  changes  in  solar  illumination  induce  abrupt  changes  in  ground  temperature  (>  IK  in  the 
first  minute),  and  the  cooling  or  warming  continues  at  a  slower  rate  for  over  an  hour  with  a 
total  change  in  temperature  of  3  to  4K. 


ca 
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Figure  18.  Diurnal  Variation  of  Surface  Air,  Ground,  Surface  Soil,  and  Barley 
Temperature  Under  100%  Cover  by  Stratus  Clouds,  0.33-1.00  km. 
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Figure  19.  Diurnal  Variation  of  Temperatures  Under  10%  Cover  by  Stratus  Clouds. 
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Figure  20.  Diurnal  Variation  of  Temperatures  Under  40%  Cover  by  Stratus  Clouds 
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TEMPERATURE  DIFFERENCES  (K)  y  TEMPERATURE  (K) 


TIME  (HRS  AFTER  MIDNIGHT) 

Diurnal  Variation  of  Temperatures  Under  70%  Cover  by  Stratus  Clouds. 


TIME  (HRS  AFTER  MIDNIGHT) 

Figure  22.  Differences  in  Ground  Temperatures  Between  Stratus  Clouds  (0.33-1.00  km) 
and  Cloud-Free  Conditions. 
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The  effect  of  the  large  contribution  from  skyshine  during  the  day  is  seen  in  the 

temperature  difference  curves  of  Figure  22.  The  ground  temperature  curves  under  fractional 
cover  by  stratus  clouds  do  not  fall  between  the  cloud-free  and  100%  cloud  cover  limits, 

especially  during  midday.  The  reason  is  as  follows:  The  cloud  lattice  always  has  a  cloud 

directly  overhead  which  occults  the  sun  between  9  a.m.  and  3  p.m..  Figure  6.  Thus  these 

clouded  skies  are  always  cooler  than  the  cloud-free  sky.  In  addition,  the  skyshine  contribution 
increases  with  increasing  cloud  cover.  Since  the  10%  cloud  cover  curve  has  the  least  warming 
due  to  skyshine  and  still  has  no  midday  solar  illumination,  the  ground  it  shadows  is  coolest. 


* 
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4.  RECOMMENDATIONS 


This  study  has  demonstrated  the  necessity  of  including  cloud  shadowing  in  the  SWOE 
scene  model.  It  has  also  produced  a  series  of  modules  which  can  serve  as  a  basis  for 
development  of  *»  comprehensive,  dynamic  cloud  shadowing  model.  Model  upgrades  and 
validation  are  needed  to  meet  SWOE  performance  requirements.  SSI  recommends  the  following 
upgrades  be  implemented: 


•  Introduce  spectral  emissivities  and  albedos  into  the  atmospheric  calculations 
and  the  energy  budget  model. 

•  Use  the  SSI  reflectance  model(^)  to  account  for  the  surface  scattering  angular 
dependence,  especially  with  regard  to  the  absorption  of  solar  radiation. 

•  Generalize  the  cloud  fields  to  allow  translation  (wind),  growth  and 

inhomogeneity  of  clouds.  Previously  developed  cloud  distribution  modules,  such 
as  those  developed  by  TASC(^)  and  ONTAR(^),  could  be  introduced  into 
the  shadowing  model. 

•  Calculate  skyshine  and  solar  irradiance  using  the  SSI  cloud  radiance  modelX11) 

•  Account  for  the  finite  size  of  the  sun. 

•  Add  Deardorffs  water  conservation  equations^)  to  the  terrain  model. 

•  Replace  the  simple  single-layer  vegetation  model  currently  in  TERRAIN  with 
Sparta’s  forest  model. 
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